UNCLASSIFIED 


AD  NUMBER 


AD843586 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  Sep  1968. 
Other  requests  shall  be  referred  to  Office 
of  Naval  Research,  Washington,  DC  20360. 


AUTHORITY 


ONR  Notice,  27  Jul  1971 


THIS  PAGE  IS  UNCLASSIFIED 


A® 843586 


HPC  68-68 


i 

MONSANTO/WASHINGTON  UNIVERSITY 
ONR/ARPA  ASSOCIATION 

CRACK  GROWTH  IN  UNIDIRECTIONAL  COMPOSITES 
UNDER  REPEATED  LOADINGS 

Edward  M.  Wu 


September  1968 


PROGRAM  MANAGER 
ROLF  BUCHDAHL 


This  document  is  subject  to  special  export  controls  end  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be  made  only  with  prior  approval  of  the  Director 
of  Material  Sciences,  Office  of  Navel  Research.  LOASH  o20  3(c>O 

MONSANTO  RESEARCH  CORPORATION 

A  SUHIOIAIT  OF  MONSANTO  COMtANT 

too  N.  LINDBERGH  BOULEVARD  ST.  LOUIS.  MISSOURI  ASIM 


D  D  C 

iTPfp.fi3  pn 

NOV  22 1968 

3tbiJ  I?  ’■ 

V 


r 


•  •V 


Hi 


w*n  ttnr 


WTIOI  □ 


WT  KSTIfllj 


NOTICES 

When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Govern* 
ment  procurment  operation,  the  United  States  Government  thereby 
incurs  no  responsibly  nor  any  obligation  whatsoever;  and  the  fact  that 
the  Government  may  have  formulated,  furnished,  or  in  any  way  supplied 
the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded 
by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 

DDC  release  to  CFSTI  is  not  authorized. 


AD-  71* 


HPC  68-63 


♦  , 


J 

/ 


MONSANTO/WASHINGTON  UNIVERSITY  ASSOCIATION 
Sponsored  by  ONR  and  ARPA 

Development  of  High  Performance  Composites 

CRACK  GROWTH  IN  UNIDIRECTIONAL  COMPOSITES 
UNDER  REPEATED  LOADINGS 

Edward  M.  Wu 
September  1968 

Rolf  Buchdahl,  Program  Manager 


Monsanto  Research  Corporation 
800  North  Lindbergh  Boulevard 
St.  Louie ,  Missouri  63166 


This  paper  was  prepared  under  the  sponsorship  of  the  Advanced 
Research  Projects  Agency,  Department  of  Defense,  through  a 
contract  with  the  Office  of  Naval  Research,  N00014-67-C-0218 
(formerly  N00014-66-C-0045) ,  ARPA  Order  No.  873,  ONR  Contract 
Authority  NR- 356-484/4-13-66 ,  "Development  of  High  Performance 
Composites . " 


FOREWORD 


The  research  reported  herein  was  conducted  by  the  staff  of 
the  Monsanto/ Washington  University  Association  under  the  sponsor 
ship  of  the  Advanced  Research  Projects  Agency,  Department  of 
Defense,  through  a  contract  with  the  Office  of  Naval  Research, 
N00014-67-C-0218  (formerly  N00014-66-C-0045),  ARPA  Order 
No.  873,  ONR  contract  authority  NR  356-484/4-13-66,  entitled 
"Development  of  High  Performance  Compositea.  " 

The  prime  contractor  is  Monsanto  Research  Corporation. 

The  Program  Manager  is  Dr.  Rolf  Buchdahl  (phone  314-694-4721). 
The  contract  is  funded  for  $5,000,000  and  expires  30  April 


1970. 


CRACK  GROWTH  IN  UNIDIRECTIONAL  COMPOSITES 
UNDER  REPEATED  LOADINGS 

Edward  M.  Wu* 

ABSTRACT 

A  phenomenological  description  of  the  slow  growth 
of  longitudinally  oriented  cracks  in  an  orthotropic 
plate  is  explored.  Detail  microscopic  mode  of  crack 
extension  is  dominated  by  crack  skipping  around  the  fibers. 
This  crack  skipping  contributes  to  the  increase  of  re¬ 
sistance  to  subsequent  crack  growth.  Mathematical 
models  are  used  to  relate  experimental  results  of  crack 
growth  under  repeated  loadings. 


^Materials  Research  Laboratory,  Washington  University, 
St.  Louis,  Mo. 


Introduction 


The  rupture  process  of  structural  materiel  due  to  repeated  loadings  (fatigue) 

Is  generally  classified  into  three  stages.  These  stages  are:  nucleatton  of  flaws,  slow 
growth  of  flaws  to  cracks  of  Griffith  dimension,  and  finally,  rapid  propagation  of 
cracks  to  catastrophic  failure.  It  is  the  general  consensus  that  a  majority  of  the  life 
time  of  a  structure  is  consumed  In  the  second  stage  of  quasi-equilibrium  slow  growth 
of  cracks  [1].  Henee,  an  understanding  and  characterization  of  the  slow  growth  of 
cracks  under  repeated  loading  will  contribute  to  a  more  reliable  prediction  of  the 
lift  of  a  structure.  Previous  Investigators  have  attempted  to  characterize  this  crack 
growth  stage  by  employing  concepts  equivalent  to  limiting  stress  (Frost  &  Dugdale  [2], 
Liu  C31),  limiting  strain  (Head  [4]),  or  energy  stability  (Paris  [5],  Uu  [6],  Cotterell 
[7]).  These  Investigations  are  limited  to  statistically  homogeneous  isotropic  mate¬ 
rials.  In  view  of  the  recent  Increase  of  composite  materials  in  structural  application, 

It  is  desirable  to  extend  these  approaches  to  anisotropic  materials.  As  an  initial 
attempt,  the  slow  growth  of  cracks  in  orthotropic  plates  is  explored.  It  is  hoped  that 
this  work  will  provide  a  foundation  toward  an  understanding  of  general  classes  of  aniso¬ 
tropic  and  laminated  bodies. 

Among  the  aforementioned  approaches  in  characterizing  slow  crack  growth 
In  Isotropic  materials,  the  energy  stability  or  crack  driving  force  concept  appeared 
to  provide  the  most  sjtlsfectory  description  over  a  broad  spectrum  of  crack  growth 


velocities.  Furthermore,  this  approach  also  provides  a  link  between  crack  growth  by 
repeated  loadings  and  the  Griffith  theory  of  fracture.  This  general  approach  is  adopted 
for  our  analysis  of  crack  growth  in  orthotropic  plates. 


Phenomenological  Description  of  Crack  Growth 


In  this  investigation,  a  longitudinally  oriented  crack  in  a  unidirectional  fiber- 
reinforced  composite  is  examined.  An  attempt  is  made  to  relate  the  slow  growth  in 
static  fracture  to  that  under  repeated  loadings.  The  Griffith-type  condition,  generally 
employed  in  linear  fracture  mechanics,  essentially  assumes  that  the  onset  of  cata¬ 
strophic  crack  propagation  is  a  function  of  critical  energy  (or  driving  forces  as  inter¬ 
preted  by  Irwin  [8] ).  Time  effect  is  ignored  and  crack  geometry  is  assumed  to  be 
independent  of  the  magnitude  and  history  of  loading.  In  the  fracture  of  composites, 
these  assumptions  must  be  examined  critically  and  modified  whenever  necessary. 


The  propagation  pattern  of  a  longitudinally  oriented  crack  in  a  unidirectional 
fiber-reinforced  composite  differs  significantly  from  that  In  a  homogeneous  isotropic 
medium  [9].  The  gross  propagation  trajectory  In  a  unidirectional  composite  is  eo- 
linear  with  the  original  crack  both  under  symmetric  and  skew-symmetric  loadings. 
However,  the  latter  has  not  been  observed  in  isotropic  materials.  In  fact,  these  pat¬ 
terns  uniquely  coincide  with  the  opening  and  forward  shear  modes  envisioned  by  Irwin. 
Treating  a  material  as  statistically  homogeneous  and  orthotropic,  its  strain  energy  re¬ 
lease  rates  or  the  equivalent  parameters— stress  Intensity  factors— can  be  easily  evalu¬ 
ated.  The  critical  stress  intensity  factors  can  be  determined  experimentally  and  they 
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proved  to  form  a  vary  satisfactory  crltarla  [9].  To  axtand  this  treatment  to  crack 
growth  under  repeated  loadings,  a  closer  examination  of  crack  growth  resistance  Is 
required. 

Fracture  experiments  in  composites  [10]  indicated  that  a  significant  amount 
of  stable  crack  growth  took  place  prior  to  fast  propagation.  Inasmuch  as  the  slow 
stable  growth  and  fatigue  are  time-dependant  processes,  the  associated  anisotropic 
viscoelastic  rupture  process  must  be  characterized  and  the  experimentally  determined 
parameters  clearly  identified. 

The  effect  of  inhomogeneity  due  to  the  fibers  must  also  be  better  understood. 
Although  the  gross  propagation  trajectory  of  the  crack  is  straight,  microscope  exami¬ 
nation  revealed  that  the  crack  propagation  pattern  is  dominated  by  crack  skipping 
across  the  fibers  [11].  A  consequence  of  this  crack  skipping  phenomenon  is  that  the 
amount  of  absolute  crack  extension  increases  the  resistance  of  subsequent  crack  growth 
and  thus  introduces  a  strorg  geometric  effect  on  the  crack  growth  history. 

In  repeated  loadings,  another  Important  consideration  Is  the  local  modification 
of  material  In  the  neighborhood  of  the  crack  tip.  This  material  modification  may  take 
the  form  of  strain  hardening  in  the  case  of  crystalline  materials  and  of  cold  drawing 
in  the  case  of  amorphous  polymers.  Thus,  the  three  primary  effects  to  be  considered 
in  the  description  of  crack  propagation  under  repeated  loading  in  composites  are: 
first,  the  viscoelastic  rupture  process;  second,  the  crack  growth  history;  and  third,  the 
local  modification  of  the  material  near  the  crack  tip.  The  viscoelastic  rupture  of 
Isotropic  homogeneous  material  was  studied  by  Halpin  [12].  The  geometric  effect 
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of  crack  growth  history  Is  traatad  In  data  I  i  here,  whila  tha  remaining  two  effects  for 
composites  ara  laft  for  future  invasttgatlon. 

Analysis  of  tha  Geometric  Effect  of  Crack  Growth 

In  this  section/  the  geometric  affect  of  crack  skipping  cumulated  In  the  course 
of  slow  growth  is  first  examined.  Its  role  In  relation  to  the  rata  of  crack  growth  undar 
repeated  loading  is  than  explored. 


Tha  crack  skipping  phenomena  Is  attributable  to  tha  strong  orthotropy  In 
strength  and  stiffness  of  tha  composite  [10].  In  Isotropic  materials,  subcrltical  flaws 
can  be  envisioned  to  distribute  randomly  over  the  body.  Under  the  Intensified  stress 
field  In  the  neighborhood  of  the  main  crack  tip,  the  flaw  which  Is  normal  to  the  maxi¬ 
mum  tensile  stress  (cr^  Fig.  la)  will  extend.  The  main  crack  then  propagates  by 
joining  with  flaws  In  this  direction.  Under  many  loading  conditions,  a.g.,  skew- 
symmetric  loads,  this  direction  of  propagation  Is  not  co-l Inear  with  the  crack.  In 
unl directionally  reinforced  composites,  however,  the  subcrltical  flaws  are  randomly 
distributed  but  oriented  along  the  direction  of  the  reinforcing  fibers.  In  this  case,  the 
flaws  which  are  oriented  normal  to  the  maximum  transverse  stress  (O’  Fig.  lb) 
are  most  likely  to  extend.  Due  to  elastic  orthotropy,  the  direction  maximum  trans¬ 
verse  stress  Is  In  general  not  co-l  Inear  with  the  crack.  Ergo,  when  these  flaws  are 
joined  with  the  main  crack,  It  then  gives  the  appearance  of  the  crack  skipping,  while 
leaving  strings  of  fiber  in  the  path  of  propagation.  On  the  other  hand,  because  the 
strength  of  the  fiber  Is  usually  much  higher  than  the  matrix,  the  crack  is  channeled 
in  the  general  direction  of  the  matrix.  According  to  the  elastic  stress  analysts,  this 
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crack  skipping  should  ba  random  when  subjected  to  pure  tension  (symmetric  load) 
and  It  should  be  oriented  In  the  direction  of  shear  when  subjected  to  pure  shear  (skew- 
symmetric  load)  [10]. 

Figs.  2  are  closed-up  photographs  of  the  crack  skipping  modes.  Fig.  2a  shows 
a  crack  which  skipped  randomly  under  tension.  Fig.  2b  shows  a  crack  which  skipped 
in  the  direction  of  the  shear  force  and  Fig.  2c  shows  a  crack  which  reversed  Its  direc¬ 
tion  of  skipping  whan  the  shear  force  was  reversed.  Figs.  2d  and  2e  are  further  evi¬ 
dences  of  crack  skipping.  Fig.  2d  shows  a  specimen  broken  in  tension  with  randomly 
oriented  strings  of  fiber  still  attached  to  both  sides.  Fig.  2e  shows  a  specimen  broken 
under  a  combined  tension  and  shear  now  with  strings  of  fibers  oriented  in  the  direction 
of  shear  force. 

To  assess  the  effect  of  these  strings  formed  by  unbroken  fibers,  one  can  con¬ 
ceptually  remove  these  stringers  and  replace  them  by  singular  forces  [10].  In  the 
course  of  summing  up  the  effect  of  these  forces,  several  difficulties  arise.  First  of  all, 
the  number  of  stringers  formed  cannot  be  conveniently  evaluated  experimentally. 
Secondly,  the  total  effect  of  the  stringers  cannot  be  obtained  by  simply  using  the  singu¬ 
lar  force  solution  as  Green's  function  integrating  over  the  length  of  crack  extension 
because  a  singular  force  exerted  at  the  crack  tip  is  also  involved.  Both  of  these  prob¬ 
lems  can  be  alleviated  by  assuming  that  the  effect  of  the  stringers  Is  equivalent  to 
normal  and  shear  forces  distributed  over  the  length  of  crack  growth.  These  distributed 
forces  are  functions  of  the  stiffness  of  the  reinforcing  fibers  as  well  as  the  relative  dis¬ 
placement  between  the  top  and  bottom  of  the  crack  surfaces.  These  forces  can  be 
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evaluated  if  the  precise  geometry  of  crack  skipping  is  known.  However,  since  these 
can,  at  best,  be  predicted  in  a  statistical  manner,  it  would  be  far  more  practical  to 
assume  that  they  are  uniform  as  illustrated  schematically  in  Fig.  3.  Their  effective 
magnitude  p  and  t  shall  be  determined  experimentally.  The  effect  of  these 
stringer  forces  are  analyzed  first  and  their  combined  effect  with  externally  applied 
forces  shall  be  arrived  at  by  superposition.  The  linear  elastic  boundary  value  problem, 
shown  in  Fig.  4,  can  be  solved  by  the  method  described  in  Ref.  [131.  Consider  the 
shear  force  alone  first  in  the  mapped  plane  (Fig.  4b),  the  boundary  conditions  are:* 

4  -c 


4 

't'si 

•»  c 

.  —'i 

•'*\ 

■A 

-  l: 

%  r. 

where  or  is  a  point  on  the  unit  circle  and  0|  Is  arbitrarily  chosen  as  the  starting 
point. 


Substituting  Eq.  (1)  into  Eq.  (30)  in  Ref.  [13]  and  upon  evaluating  the  inte- 
g  als,  the  following  result  can  be  arrived  at: 

=  *-  ~~J(i  - h)\ U'\'  — ■»  f-  -  iQ^I— ~ - 1 — V 

4  V  \  Vr/-C)  |>  cr(*-Cu 


-  r1 

-  W  \'2 


*The  same  notations  as  in  Ref.  [1  3]  are  used  here. 
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Proceeding  in  a  similar  manner,  the  results,  for  uniformly  distributed  no -mat  forces 
(Fig.  3b),  are: 


r>a  f, 


-v  iOAo 


_L\r  ^  —  u-iii ,,  ,j.  _l —  - .  ^ .  1 

rrVr^V  ltT  '  V-T.V'  v;i 

-  -v*  r1 

l 

(3) 


v,V) 


V  A  •  l 


Using  these  complex  potentials  ) ,  the  stress  distribution  can  be  deter¬ 

mined  for  standard  relations  given  by  Lekhnitsklf  [14). 


Since  both  flaw  nucleatlon  and  crack  growth  initiate  from  the  vicinity  of 
the  crack  tip,  the  solution  can  be  specialized  for  this  region  by  taking  the  limit  of 
when  ^  approaches  unity  (  ^  *  1  being  the  crack  tip).  Performing  the 

<) 

limiting  process  on  Eq.  (2)  and  (3),  and  noting  that 


the  following  results  can  be  obtained,  for  a  distributed  shear, 


t ,  over  A  a: 


b,«\  •  X  iu _ Q{.  -*/  -  *>.'  \ 

and  for  a  distributed  tension,  p  ,  over  a,  a 


(4) 


lv,  t  .  cc: 

1  s  <v  )  <s) 

From  these  results,  together  with  the  relationship  given  by  Eq.  (3)  in  Ref.  [8],  the 
stress  intensity  factors  can  be  obtained  directly.  Using  the  subscript  s  to  denote  the 
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stress  intensity  factors  associated  with  the  stringers,  for  the  distributed  tension  and 
shear,  respectively,  they  are: 


k  '' 
-- 


££  o-  " c.  ( 

Ti  k  rv  > 


2  C 


(\  ‘  cr  -  V 

ii  '  tv  ) 


(6) 


As  stated  earlier,  the  macroscopic  crack  propagation  modes  are  precisely 
mode  1  and  mode  II.  Consequently,  the  strain  energy  release  rates  can  be  directly 
related  to  the  stress  intensity  factors  by  (Cf.  Ref.  [15] ): 

X  -  K  [*){*  ?£  * 


Thus,  for  the  cases  under  consideration,  the  strain  energy  release  ^j.  and 

2 

k.  are  equivalent  parameters.  For  the  sake  of  simplicity  in  computation,  stress 
intensity  factors  are  used  in  this  study. 


If  kj  are  used  to  denote  stress  intensity  factors  when  the  effect  of  stringers  is 
ignored,  then  the  actual  stress  intensity  factors  at  the  crack  tip  can  be  obtained  by 
superposition  in  the  form: 


k  =  k  +  k 

1  1  IS 


k  =  k  *4"  k 
*2  2  k2s 


(8) 


Sinca  both  k^  and  k2>  ara  nagativa  (Eq.  6),  tha  stringers  formad  by  crack 
iklpplng  tend  to  retard  tha  crack  growth,  or  tantamountly,  they  Increaie  tha  material 
resistance  to  crack  growth.  From  Eq.  (6)  and  (8)  it  can  be  seen  that  the  longer  the 
crack  growth  {  A  a),  the  greater  tha  resistance  to  crack  growth.  Insofar  as  crack 
propagation  by  repeated  loading  is  a  quasi -equilibrium  crack  growth  process,  the 
growth  resistance  in  such  composite  should  increase  with  the  loading  cycles.  This 
phenomena  can  be  characterized  by  modifying  the  Cotterell-Kraffs  [6,  16]  model 
for  fatigue  crack  growth  in  isotropic  materials. 

In  Kraffs  model  [16],  the  Instability  of  a  crack  is  defined  when  the  slope  of 
the  crack  resistance  denoted  as  R  is  tangent  to  the  driving  force  denoted  as  . 

2 

Since  ^  and  (k)  differ  by  a  constant  which  Is  a  function  of  the  material 

(Eq.  7),  the  same  relationship  exists  as  Indicated  schematically  in  Fig.  5.  In  the  first 

2 

cycle  of  loading,  the  critical  state  is  k  .In  the  second  cycle,  due  to  crack  growth 

C1 

■A  the  growth  resistance  increases  as  indicated  as  R2  which  raises  the  critical 

2 

point  to  k  .  In  a  similar  manner,  the  cumulated  crack  growth  Aan  raises  the 
c2  -  2 

2 

critical  state  to  a  still  higher  point  k  .  This  phenomena  differs  from  those  observed 

c3 

in  the  isotropic  case  in  that  the  experimentally  observed  crack  resistance  curve  in 
isotropic  materials  and  the  critical  stress  intensity  factor,  for  all  practical  purposes, 
are  Independent  of  crack  length  and  growth  history.  Although  in  metals  the  crack 
resistance  curve  changes  somewhat  in  the  initial  cycling  drive  to  strain  hardening,  it 
stabilizes  in  subsequent  cycles. 
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In  order  to  modify  Cotterell's  [6]  model  for  fatigue  crack  propagation  for 
this  class  of  composites,  It  Is  necessary  to  hypothesize  thet  the  shape  of  the  crack 
resistance  curve  R  Is  unique  and  that  It  Is  Independent  of  crack  length  but  Is  a 
function  of  the  mode  of  the  externally  applied  load  P(0)*,  history  of  crack  growth 
A  a,  time  t,  and  temperature  T,  i.e., 

r.  "  U  m  (9) 

Projecting  the  shape  of  the  crack  resistance  curve  R  schematically  In 
Fig.  6  for  one  mode  of  loading  (8  =  constant)  a  crack  propagation  model  for  repeated 
load  can  be  obtained.  For  a  crack  of  Initial  length  aQ  under  a  repeated  load  of  con¬ 
stant  amplitude,  so  the  crack  growth  would  follow  the  appropriate  crack  resistance 
curve  to  length  a^  until  the  load  Is  released.  In  the  second  cycle,  the  crack  would 
follow  a  higher  resistance  curve  (higher  because  of  the  growth  (a.  -  a>)  and  grows 
to  a2>  Similarly,  In  the  third  cycle  the  crack  would  follow  a  still  higher  resistance 
curve  and  grows  to  a^.  This  process  is  continued  until  In  the  last  cycle  Is  tangent  to 
the  load  curve  at  which  time  rapid  fracture  occurs.  Accordingly,  If  the  function  sug¬ 
gested  by  Eq.  (9)  Is  known,  the  number  of  cycles  to  fracture  can  be  estimated. 

Since  the  time-temperature  effect  Is  kept  constant  in  this  study,  only  the  first 
two  parameters  in  Eq.  (9)  are  considered.  Since  crack  growth  resistance  is  a  material 
property,  the  hypothesis  that  the  crack  resistance  curve  is  unique  has  to  be  tested 
experimentally.  If  the  resistance  curve  Is  Indeed  unique,  then  the  growth  history 

*0  is  the  angle  between  the  load  vector  and  the  crack;  i.e.,  0  = 

0=0  shear. 


is  tension 
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offset  can  ba  determined  using  the  parameters  glvan  In  Eq.  (6). 

Exparlmantal  Observations 

The  goals  In  the  experimental  observation  are:  first,  to  Investigate  the 
Increase  In  the  critical  point  of  instability  both  under  tension  and  under  shear 
(0  =  ^  and  0=0,  respectively,  under  the  notation  used  in  Eq.  9);  second,  to 
explore  whether  the  crack  rasistanca  curve  Is  unique  whan  the  mode  of  load  Is  the 
same  (9  *  constant);  third,  to  check  whether  the  mathematical  models  are  ade¬ 
quate  In  characterizing  the  crack  growth  history;  and  finally,  to  estimate  the  effec¬ 
tive  tension  and  shear  forces  (p  and  t  In  Eq.  6)  Induced  by  the  stringers. 

The  materials  used  in  the  experiments  are  3M  Scotch-ply  Type  1002  uni¬ 
directional  fiberglass  reinforced  plastic  plates  0.05  inch  thick.  In  each  specimen, 
one  straight  crack  Is  cut  parallel  to  the  reinforcing  fibers.  Rectangular  specimens 
6  inches  wide  by  10  inches  long  were  used  for  symmetric  loading  tests  in  which  tension 
Is  applied  perpendicular  to  crack  and  fiber.  Cantilever  beam  specimens  4  inches  high 
by  10  Inches  long  were  used  for  skew-symmetric  loading  tests.  It  was  astablished  in 
Ref.  [14]  that  the  cantilever  beam  was  equivalent  to  a  pure  shear  fracture  test  when 
a  crack  was  located  along  the  neutral  axis.  The  experimental  set-up  and  procedure 
are  similar  to  those  desceibed  in  Ref.  [9].  Briefly,  the  procedure  consists  of  using  a 
motor-driven  camera  to  record  the  crack  length  in  the  process  of  loading.  An  event 
marker  In  the  load  recording  chart  Is  synchronized  with  the  load  to  the  crack  length 
recorder  on  the  film.  From  these  records,  point  wise  crack  growth  resistance  can  be 
computed  from  which  a  continuous  curve  can  be  drawn.  The  specimens  were  subjected 
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to  approximately  four  cyelas  ai  ihown  schematically  in  Pig.  5.  In  the  first  cycle, 
load  was  increased  while  the  crack  growth  was  closely  monitored.  When  the  com¬ 
bination  of  load  and  oack  propagation  appeared  to  have  reached  or  slightly  passed 
the  instability  point  (the  point  of  tangency  of  the  driving  force  and  the  resistance 
curve),  the  specimen  was  rapidly  unloaded.  This  was  repeated  for  the  second  and 
third  cycles  and  in  the  fourth  cycle  the  specimen  was  loaded  all  the  way  to  fracture. 

In  the  symmetric  load  tests,  the  loading  program  was  zero  to  maximum  tension  (up 
to  instability)  for  four  cycles.  In  the  skew-symmetric  load  tests  the  loading  program 
was  zero  to  maximum  shear  In  alternate  directions  for  four  cycles.  The  typical  crack 
resistance  curves  thus  obtained  for  tension  and  shear  respectively  are  shown  In  Figs.  7 
and  8.  The  stress  intensity  factors  kj  are  computed  where  the  stringer  forces  are 
not  taken  into  consideration.  They  are  computed  by  the  relationships  [8,  16] 
k j  *  /a  and  k^  =  jjj-  /a  where  <r  is  the  tension,  a  is  the  half  crack  length 

and  Q  the  shear  force  A  the  transverse  cross  section  area. 

The  critical  values  plus  other  pertinent  Information  is  tabulated  in  tables  I 
and  II  for  tension  and  shear, respectively.  From  Figs.  7  and  8  it  can  be  noted  that, 
with  experience,  the  specimens  were  loaded  very  close  to  or  sometimes  past  the  critical 
point  In  the  repeated  loading  process.  This  exparienee,  however,  was  not  gained  with¬ 
out  cost;  the  missing  specimen  numbers  in  tables  I  and  II  were  broken  in  this  learning 
process. 


Both  of  the  figures  and  the  tabulated  values  show  clearly  that  the  critical  point 


rose  in  repeated  loadings. 
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To  explore  whether  the  creek  resistance  curve  it  unique,  the  crack  resistance 
curvet  are  normalized  and  compared.  Fig.  9  represents  the  normalized  data  of  the 
crack  resistance  data  from  four  specimens,  subjected  to  tension.  Fig.  10  represents 
the  data  from  four  specimens  subject  to  shear.  Within  the  realm  of  our  experimental 
scatter,  the  crack  growth  resistance  curves  appear  to  be  unique.  However,  significant 
differences  exist  between  the  resistance  curve  for  tension  and  shear  which  are  apparent 
when  they  are  compared  in  Fig.  11.  It  can  be  concluded  then  that  the  form  of  the 
crack  growth  resistance  curve  is  independent  of  crack  length  but  dependent  on  the 
mode  of  loading. 


To  examine  the  validity  of  the  mathematical  model  for  characterizing  crock 
skipping,  Eq.  (8)  can  be  rearranged;  together  with  Eq.  (6)  we  have, 

b,'»  V 

(10) 

u:  .  k,. 

If  the  mathematical  models  are  adequate,  the  stress  intensity  factors  k^  and 

k.  should  be  constants  and  Eq.  (10)  in  a  graph  of  k,  versus  i  '/Vcu 
should  be  straight  lines.  Oats  from  tables  I  and  II  presented  in  this  form  are  shown 
in  Figs.  12  and  13  .  Here  too,  within  our  experimental  scatter,  the  mathematical 
models  provide  very  close  characterization  indeed.  The  intercepts  in  the  ordinates 
should  be  the  critical  k^  and  k^  for  static  fracture  tests.  This  compared  very 
favorably  with  previously  obtained  results  [8].  Finally,  the  effective  tension  and 
shear,  p  and  t,  can  be  evaluated  easily  from  Figs.  12  and  13  .  The  slopes  of  the 
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2P  2t 

straight  lines  are  *  570  and  —  *  1560  respectively,  which  give  the  value* 
p  «  900  psi  and  t  *  2500  psl.  It  Is  worth  noting  that  p  Is  limited  by  the  "peel 
strength"  of  the  fiber  from  the  matrix.  The  value  evaluated  above  Is  certainly  not 
far  from  the  "peel  strength"  reported  In  other  Investigations.  Furthermore,  from  the 
same  model,  the  maximum  value  of  t  Is  limited  by  the  tensile  strength  of  the  fiber. 
In  this  case,  t  is  substantially  less  than  the  strength  of  glass  fiber  which  does  not  con¬ 
tradict  the  observation  that  the  glass  fiber  stringers  remained  Intact  after  fracture 
(Ftg.  1  and  Fig.  2). 

Conclusions 

This  Investigation  is  an  Initial  step  In  the  attempt  to  better  characterize  the 
fatigue  of  fibrous  composite  materials.  One  of  the  primary  effects,  the  geometric 
history  effect  of  crack  growth,  was  explored. 

It  was  observed  that  the  crack  resistance  curve  appeared  to  be  independent  of 
crack  length  and  of  crack  growth  history  for  a  given  mode  of  loading.  However,  the 
crack  resistance  curves  for  tension  and  for  shear  are  distinctly  different  Indicating  that 
they  are  functions  of  the  mode  of  loading. 

It  was  also  observed  that  under  small  number  of  repeated  loadings  deliberately 
designed  to  provide  large  amounts  of  crack  growth,  the  crack  growth  resistance  In- 
ereesed  with  the  number  of  cyclic  loads.  This  increase  In  crack  growth  resistance  was 
successfully  characterized  by  a  mathematical  model  in  which  the  growth  history  can 
be  related  to  static  fracture  strength.  The  magnitude  of  the  distributed  stringer  forces 
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estimated  by  means  of  this  model  also  agreed  well  with  other 
findings.  With  this  model,  the  geometrical  effect  of  crack 
growth  can  be  characterized  and  Cotterell's  mechanics  of  iso¬ 
tropic  fatigue  crack  growth  can  be  modified  for  unidirectionally 
reinforced  composites. 
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